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Численное решение краевых задач для уравнения Пуассона методом точечных 
источников поля* 


С. Ю. Князев, Е. Е. Щербакова, А. А. Енгибарян 


Целью работы является повышение эффективности одного из перспективных методов решения краевых задач 
для уравнений эллиптического типа — метода точечных источников поля. В зарубежной литературе он назы- 
вается методом фундаментальных решений. В настоящее время он используется в первую очередь для реше- 
ния уравнения Лапласа. Предложено несколько вариантов численного решения краевых задач для уравнения 
Пуассона с использованием метода точечных источников поля. Применение этого метода к решению неодно- 
родных уравнений, таких как уравнение Пуассона, приводит в большинстве случаев к резкому возрастанию 
численной погрешности, что связано с ошибками при нахождении частного решения уравнения Пуассона. 
Правая часть уравнения Пуассона аппроксимируется двумерным тригонометрическим многочленом (при ре- 
шении двумерных краевых задач), после чего становится возможным получение частного решения, необходи- 
мого для решения исходной краевой задачи методом точечных источников поля. Результаты тестирования 
предложенного способа численного решения краевых задач для уравнения Пуассона свидетельствуют о его 
эффективности, так как позволяют получать решение с относительной погрешностью 10-5 при минимальных 
затратах машинного времени. Разработанная методика численного решения краевых задач для уравнения 
Пуассона может быть использована при моделировании физических полей в технических устройствах различ- 
ного назначения. 

Ключевые слова: уравнение Пуассона, уравнения эллиптического типа, метод точечных источников поля, 
метод фундаментальных решений. 


Введение. Ряд задач, имеющих важное практическое значение, сводят к решению уравнения 
Лапласа с условиями, задаваемыми на границе с достаточно сложной конфигурацией. Одним из 
эффективных методов численного решения краевых задач для однородных уравнений эллиптиче- 
ского типа с постоянными коэффициентами, в первую очередь для уравнения Лапласа, является 
метод точечных источников поля (МТИ) [1-4]. Преимуществом МТИ является его простота и не- 
сколько меньший объём вычислений в сравнении с традиционными численными методами решения 
граничных задач, такими, например, как метод конечных элементов (МКЭ). Применение МТИ может 
быть оправдано также при решении краевых задач для неоднородных уравнений эллиптического 
типа, например, при решении уравнения Пуассона [5-9]. 


Основная часть. Рассмотрим область © < В?, в которой решается уравнение Пуассона: 


ДИ =Е(г), (1) 
с условиями на границе д, соответствующими задаче Дирихле: 
И= д. (2) 


Пусть для уравнения (1) известно частное решение (/, (г) (не обязательно с однородными 
условиями на границе): 
ЛИ, (г) = Е (г). 
Введём функцию и (г) = И (г) - Ц, (г) . Эта функция удовлетворяет однородному уравнению 


Лапласа Ди = 0 с граничными условиями 
Чо = (9 - 4 < ; (3) 





* Работа выполнена по грантам РФФИ № 13-07-00952-а и № 14-07-00705-а. 
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Для нахождения неизвестной функции и (г) можно теперь использовать МТИ, как это опи- 


сано в [3-4]. В результате искомое решение краевой задачи найдём в виде суммы 
И (г) =и(г) +, (г) . Таким образом, для успешного решения краевой задачи для уравнения Пуас- 


сона необходимо сначала найти его частное решение. При этом могут использоваться различные 
ПОДХОДЫ. 


Как известно, правая часть ! (г) в уравнении Пуассона ЛИ = # (г) пропорциональна плот- 
ности заряда р(г) в точке с координатой г: / (г) =Кр(г) . Коэффициент пропорциональности за- 


висит от используемой системы единиц и, следовательно, от конкретного представления фунда- 
ментального решения для однородного уравнения (в данном случае уравнения Лапласа). Если фун- 
даментальное решение, то есть потенциал поля, созданного в точке с координатой г единичным 


ы 1 

положительным зарядом, находящимся в точке с координатой В, имеет вид ме 
Е 

Е (г) =2пр(г). Для численного решения уравнения Пуассона область решения © разбивается на 


небольшие элементарные участки ди», каждый из которых имеет площадь бя. На этой площади 
находится заряд, приблизительно равный 
р (®,) 


2п 
где г, — координата центральной точки элементарного участка б%.. 





04, = 05, 


В результате частное решение уравнения Пуассона можно приближённо представить в виде 
суммы: 





где № — число элементарных участков, на которые разбита область решения ©. 

Таким образом, решаемая задача приводится к задаче численного решения уравнения 
Лапласа с помощью МТИ. 

Рассмотрим конкретный, тестовый пример решения краевой задачи для уравнения Пуас- 
сона. Пусть правая часть в уравнении Пуассона 


Г (г) =2(х? +у?)соз (ху) - (хЗу - ху*)т(ху). (4) 
Область решения представляет собой квадрат, сторона которого равна единице: 
@={-0,5<х<0,5, -0,5<у<0,5}. (5) 


Решается задача Дирихле с граничными условиями, соответствующими точному решению 
задачи 


И (г) = ху ят(ху). (6) 

Условия (4) — (6) соответствуют условиям задачи, решаемой в [10] путём разложения по 
собственным функциям уравнения Гельмгольца. При численном решении уравнения Пуассона с 
условиями (4) — (6) с помощью МТИ заряды, моделирующие искомое поле для уравнения Лапласа 
располагались на окружности радиусом, равном 1,25. Заряды бд, внутри области © располагались 


в узлах равномерной прямоугольной сетки, покрывающей область ©. В узлах этой сетки на границе 
области © располагались узлы коллокации, число которых соответствовало числу зарядов для урав- 
нения Лапласа. 

На рис. представлен график зависимости погрешности численного решения тестовой задачи 
от полного числа зарядов, моделирующих искомое поле (тонкая сплошная линия 1). 
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Рис. Зависимость погрешности численного решения от количества зарядов и узлов в прямоугольной (линии 1; 2) и 
круговой (линия 3) области, полученные различными способами: моделированием частного решения зарядами внутри 
области © (линия 1) и использованием разложения правой части уравнения Пуассона в двумерный тригонометрический 
ряд (линии 2; 3) 


Анализ данных, представленных на рис., показывает, что величина погрешности МТИ об- 
ратно пропорциональна числу зарядов (точнее, = ^^ №"). Это свидетельствует о достаточно мед- 
ленном убывании погрешности МТИ с ростом числа зарядов. В то же время очевидно, что уже при 
небольшом числе зарядов, равном 25, погрешность = ^^ 103. Следовательно, в тех случаях, когда 
требования к точности решения не являются слишком строгими, применение МТИ для получения 
приближённого решения уравнения Пуассона описанным выше способом является вполне допусти- 
мым и разумным. 

Повысить точность численного решения краевой задачи для уравнения Пуассона можно, 
аппроксимируя правую часть уравнения (1) тригонометрическим многочленом. 

Известно [11], что функцию ф(х), определённую в узлах ю = хи, №, ..., хина отрезке / можно 
аппроксимировать тригонометрическим многочленом, имеющим вид 


ч(х)=>.Ф (х), 


и ви [п 5%) 


КАК = [п Х; г Хх | 


При использовании равноотстоящих узлов функции Р; (х) могут быть преобразованы втри- 








где ф, =Ф(х,), Р(х)= 


гонометрические многочлены [12]. В результате получаем 
1х м. 
Ф(х)= $ Ут (Вх +с, )} 
1= К=1 


где 








2п п 1 [2(1-1) 2х 1 
=—|А-—|, 1<А<П: С, ь - 
Ь, "| "| Пу СЕ=П 5 | р Г }(‹ з) ‚1<К<П (7) 


Аналогичным образом функцию двух переменных Е(х, У) (правая часть уравнения (1)), 


заданную в прямоугольнике со сторонами /, уи определённую в узлах прямоугольной равномерной 
сетки (х„, у " ‚1</<п,, 1< ] < п, , можно представить в виде тригонометрического многочлена 








Е(х,у) = 1 > > [5 р, п, У эп (ех + сл, ) т (Ву «вр 
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где Г, = И а параметры Б,,,Б, ,‚с»,сь, и другие определяются по формулам, аналогич- 
ным (7). 
Так как частным решением уравнения Пуассона при 
Г(х,у) = эт(Ь*х + с^ эт (Ь’у + с”) 
эп (Ь”х а с” эт (Бу ь С" 
(5%) + (5) 


ное решение уравнения (1) запишется в виде: 








является функция и (г) ‚ то, учитывая линейность уравнения, част- 


-— 
з 
> 
з 
< 
з 


о эп(Вих си, эту + с», 


лю 
— й х \2 2 
ху 1= 1 АХ =1 Ку=1 (№) + (%,) 








(8) 


На рис. представлен график зависимости погрешности численного решения тестовой задачи 
от полного числа зарядов, моделирующих искомое поле (жирная сплошная линия 2), полученный 
с использованием частного решения (8). Видно, что погрешность решения в этом случае почти на 
два порядка ниже, чем в том случае, когда поле, соответствующее частному решению, моделиру- 
ется зарядами, расположенными внутри области ©. 

При решении краевой задачи для области ©, форма которой отличается от прямоугольной, 
строится прямоугольник минимальных размеров (для минимизации численной погрешности), охва- 


тывающий область ©; функция Ё (х, у) экстраполируется за пределы области © вплоть до границ 


построенного прямоугольника. Далее находится численное решение краевой задачи описанным 
выше способом. При этом может оказаться, особенно при решении задач Дирихле, что погрешность 
численного решения для области © окажется ниже, чем для прямоугольной области. Это объясня- 
ется тем, что в области © отсутствуют участки прямоугольника, примыкающие к его вершинам, в 


которых погрешность аппроксимации функции Г (х,у) наибольшая. Эту особенность описывае- 


мого здесь метода иллюстрирует пунктирная линия 3 на рис., построенная при решении тестовой 
задачи для круговой области радиусом А = 0,5. Очевидно, что точность решения для круговой 
области почти на порядок выше, чем для прямоугольной области. При использовании зарядов, 
расположенных внутри области решения непрямоугольной формы (линия 1 на рис.), точность ре- 
шения значительно ухудшается. 

Заключение. Полученные результаты свидетельствуют о возможности эффективного использо- 
вания МТИ не только при решении краевых задач для уравнения Лапласа, но и для неоднородных 
уравнений эллиптического типа, например для уравнения Пуассона. 

Следует отметить, что предложенный алгоритм получения частного решения уравнения 
Пуассона может применяться и к другим неоднородным уравнениям эллиптического типа, например 
к уравнению Гельмгольца. Возможно также обобщение предложенного метода решения уравнения 
Пуассона для трёхмерного случая. 
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МИМЕВТСАЕ 5ОГОТТОМ ТО ВОУМОАКУ РВОВЕЕМ$ РОВ РОТ$$ОМ ЕОЦЧАТТОМ ВУ РОТМТ- 
ЗОЧВСЕ МЕТНОО* 


$. У. Кпуагем, Е. Е. ЗИспегЬаКома, А. А. Уепдфагуап 


Тре айт оГ {#65 рарег /5 пе етасепсу 1итргоуетепЕ оЁ опе о! {Те тозЕ а мапсеа {есптдиез оЁ зоМпд {те е!Ирис 
Боипаагу уаие ргоетз — {те Пе ротЕ-зоигсе тетоа ае$'дпаёеа аз {Те пдатега! зо/ивоп {есптдие т {пе 
Югеоп егаёиге. Мои! # 15 изеа рита Юг о/мтд [ар/асе едиайоп. 5еуега/ айетае питегса! 5о/ийоп$ {0 {Пе 
Боипдагу чаше ргоетз юг Ро5зоп едиайоп итд {Ве Пе/а ро/пЕ-зоигсе теЁпоа аге ргорозеа. Т5 тепоа арр/- 
савоп ю ше поппотодепеоиз едиайоп зо/ивоп, зисв аз Рой55оп едиайоп, п тозЁ сазез [еа$ ю те агатанс 
псгеазе о! {те питегса! еггог дие ю тиакКез т Розоп едиайоп зресс зо/ивоп. Тпе парЕ тетвЕГ о!Ё Ро/550п 
едиайоп 1/5 арргохтаеа Бу а Вио-@тепзюпа! {"1допотейтс роупотй! (т {те 5о/иНоп оЕЁио-@тепзюпа! Боипаагу 
уаие ргоБетз), {еп # Бесотез роз {0 оМат {Те зресйс 5о/ийоп песез5агу Юг 5оМтд ап та! Боипаагу 
уаие ргоМет Бу пе пе/а ротЕзоигсе тепод. Тпе {езНтд гези!5 о! {те ргорозед {есптдие (тру 5 ейсиепсу, ах 
теу а!ои/ оматта пе зоивоп ий а ге@вуе етог о! 10-6 аЁ тпйтит тасйте #те зрепатпа. Тйе 4еуеореа 
{есптдие о {те питегса! зо[иноп ю {те Боипаагу уа/ие ргоБетз юг Ро/55оп едиайоп сап Бе изед юг тодЕета 
рпуза! Ве/а$ т {пе епдтеейпд 4емсез о! иапоиз арр/саНопз. 

Кеуигога&: РоГ5оп едиайоп, еШрёс Боипаагу уаие ргоМетз, пе ропЁ-зоигсе те под, тепоа о! Гипаатегиа! 
50ИНОП$5. 





* Тне гезеагсН 15 допе м/й {пе зиррой: гот ВЕЕТ (дгап$ поз. 13-07-00952-а, 14-07-00705-а). 
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